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We develop a numerical algorithm for identifying approximately conserved quantities in models 
perturbed away from integrability. In the long-time regime, these quantities fully determine corre¬ 
lation functions of local observables. Applying the algorithm to the perturbed XXZ model we find 
that the main effect of perturbation consists in expanding the support of conserved quantities. This 
expansion follows quadratic dependence on the strength of perturbation. The latter result together 
with correlation functions of conserved quantities obtained from the memory function analysis con¬ 
firm feasibility of the perturbation theory. 
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I. INTRODUCTION 

A considerable interest has recently been attracted by 
integrable quantum models which, in contrast to generic 
systems, have macroscopic number of local conserved- 
quantities (CQ). Due to their presence, the isolated in¬ 
tegrable systems don’t thermalizeii^ but instead relax 
towards non-thermal steady statesii^— It has been sug¬ 
gested that such unusual steady states are fully specified 
by localiiA^ as well as quasilocal CQ. Inclusion of the 
latter ones is necessary at least in several systems which 
cannot be mapped to noninteracting particlesJ^^— The 
research on integrable systems has been motivated not 
only by such fundamental problems like that concerning 
mechanisms of thermalization/relaxation. The integrable 
systems are interesting also because they show dissipa¬ 
tionless (ballistic) energy, spin and charge transport's^— 
which might be important for the future applications. 

However, real systems are never perfect and their de¬ 
scription in terms of integrable models should be viewed, 
at best, as reasonable approximations. Therefore, it is 
important to understand the properties of systems which 
are weakly perturbed away from integrability— iSSrSS In 
the case of classical mechanics, relevant formalism has 
been developed for more than fifty years , whereas 
for quantum systems such understanding is still missing 
or, at least, remains largely incomplete. It is not evi¬ 
dent to what extend the breaking of integrability may 
be described within a single universal picture and which 
properties are specific for particular model and/or per¬ 
turbation. It is also not quite clear which hallmarks of in¬ 
tegrability disappear abruptly and which decay smoothly 
when the perturbation gradually increases. 

Recent studies allow to formulate several general 
expectations. At least at the infinite time-scale 
and for sufficiently strong perturbation the ballis¬ 
tic transport should be replaced by large but finite 
conductivity— However, after turning on the per¬ 
turbation, the system should evolve towards a quasi¬ 
steady state (prethermalizatioiJ^) which is analogous to 
the generalized Gibbs ensembles— Hence, at least for a 


finite time window [0, r] and for sufficiently weak per¬ 
turbation, the system maintains the main property of 
integrable systems, i.e., the existence of local or quasilo¬ 
cal operators which for perturbed systems are conserved 
only approximately. These quantities are not identical 
with CQ of the integrable parent model^i^i^ but rather 
they are modified by the perturbation. Finally, since 
quasilocal CQ play important role for strictly integrable 
modelsii^— they should also be included in the studies 
on perturbed systems. 

In this paper we develop an algorithm which cap¬ 
tures/verifies all these properties of the perturbed in¬ 
tegrable systems. The approach yields approximately 
conserved quantities (ACQ) which completely determine 
the long-time correlation functions of all local opera¬ 
tors supported on assumed subsystem. The algorithm 
captures cases ranging from strict integrability with lo¬ 
cal and quasilocal CQ to generic systems where CQ are 
generally nonlocal linear combinations of projections on 
eigenstates of the Hamiltonian |n)(n|. In the first part of 
this paper we show the general approach, while in the 
second part we apply it to the perturbed anisotropic 
Heisenberg (XXZ) model. We also find that the ACQ 
at weak perturbation can be described as quantities de¬ 
caying exponentially in time with the characteristic rate 
depending quadratically on the perturbation strength, 
consistent with specific findings of Ref. This result 
as well as the memory function analysis confirm non¬ 
singular behavior at weak perturbation strength and the 
feasibility of the perturbation theory for a generic class 
of integrability-breaking perturbations. 


II. GENERAL METHOD 

We study a one-dimensional tight-binding Hamilto¬ 
nian H on a lattice of L sites with periodic boundary 
conditions. We consider the space Al of local, exten¬ 
sive, translationally-invariant observables with the scalar 
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product (see Appendix A) 

{A\B) = ^{A^B) = ^Y.^*ranBmnP{Er.). (1) 

mn 

where Amn = {'^\A\n), H\n) = En\n) and the weights 
are assumed to satisfy J2nP(Bn) = 1 and p{En) > 0 
for all n. The latter assumption excludes the zero- 
temperature case but, at least in finite systems, accounts 
for the thermal states p{En) oc exp(—/3A„), where /3 is 
the inverse temperature. 

We have recently developed a procedure for identifying 
a complete set of local and quasilocal CQ in integrable 
lattice systems<^ The main step is to construct scalar 
products of all time-averaged operators {A\B), where 

1 

A = lim - / dtA{t) = Ajnn\m){n\. (2) 


which generates orthogonal set of ACQ consisting of 
Qi = ^gUsiOl. Generally, the truncation (jH) modifies 
the support of operators and transforms local operators 
Os into quasilocal ones OJ. Therefore, we split Qi into 
two orthogonal components Qi = Qf^ + such that 
= Y.si.Os\Qi)Os e Sf while = 0 for all s. 

The eigenvalues A; obtained from Eq. (| 6 ]) bear important 
information on the support of Qi (see Appendix B), 

WQiW^ = (QilQi) = Xi, = 

Carrying out the finite size scaling (L —>■ oo) of A; we dis¬ 
tinguish between local ACQ when A; = 1 for sufficiently 
large M, quasilocal ACQ when 0 < A; < 1 for any M, 
and generic nonlocal ACQ when A; —?► 0. One can also 
show (see Appendix C) that the correlation function of 
arbitrary local observables A,Bg Bff is completely de¬ 
termined by their projections on Qi, 


The identity (A|.B) = {A\B) is essential since it allows to 
distinguish between local, quasilocal and generic nonlocal 
CQ. In order to determine ACQ in the perturbed system 
one should consider a finite time-window [ 0 , r], however, 
a simple omission of the limit r —>■ oo violates the lat¬ 
ter essential relation. Therefore, we define an effective 
operator time-average with a time-scale r as 


A^ 



dt A(t) 


sin(t/T) 

TTt 


(3) 


which in spectral representation amounts to cutting off 
quickly oscillating (in time) matrix elements 

A'^ = '^e - \En - Em\) Ajnn\m){n\. (4) 

It is quite obvious that the truncated operators are ap¬ 
proximately conserved at the time-scale t t. Since 
6^(x) = 9{x), this simplified time-averaging maintains 
the property (A'^jS'^) = {A'^\B). Moreover, it becomes 
identical with the actual time-averaging over an infinite 
time-window, limT-_).oo A'^ = A, whereas for finite r it 
is related with the low-frequency spectrum of standard 
correlation functions: 


{A^\B^) = \im dw — / 

e\oJ_i 2tt L 

(5) 

Since physically interesting observables are usually 
supported on few sites only, we define a subspace 
of Al which contains operators supported on up to M 
consecutive lattice sites. The choice of interesting op¬ 
erators determines M. We introduce also the basis of 
Bff composed of operators Og which are orthonormal 
(Os|Os/) = (5s,s'. After constructing Og we solve the 
eigenproblem for the matrix ATJ^, = (OJ|OJ,) 


{A^\B^) 


{d^\Qi){Qi\E) 

I (QilQi) 

(A|Qf)(Qf|i?) 


( 8 ) 

(9) 


where the latter equation comes from identity {A\Q-j^) = 
0. Choosing A = i? as a current operator and taking 
the limit t ^ oo we recognize that Eq. dH) becomes 
the saturated Mazur boun d^^i^^ for the charge/spin stiff¬ 
ness. Note that Eq. (I 8 |) involves normalized operators 
(5i/||Q;||- Therefore, the key point is to follow how the 
supports of ACQ (and not their norms) depend on the 
time-scale r and the strength of perturbationi^ Our ap¬ 
proach gives complete set of ACQ, which are sorted from 
the most relevant local operators with A; = 1 to the least 
relevant ACQ with the smallest A;. 


III. PERTURBED ANISOTROPIC 
HEISENBERG MODEL 


Next, we apply this approach to the extended XXZ 
model 


h = jJ2 

i=i ■- 

i=i 


-iS^S^, + S-S^^ + AS]SU, 


aH' 


( 10 ) 


where are spin- 1/2 operators and the integrability 
is broken by the last term when a 7 ^ 0. We take J = 
1 as the energy unit. For concreteness, we study the 
infinite-temperature limit [p(En) = const] when Eq. ([T|) 
becomes the Hilbert-Schmidt scalar product. Then, the 
orthonormal basis of Bff is composed of operators^! 


^ UlKlg,Ugn' = SirXi, UU^ = U^U = 1 , ( 6 ) 






J + 1 
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^j+m-1-: 


m = 2,...,M (11) 
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where cr| = = ^/2Sf,a° = l,s = (si,...,Sm), 

Sj e {+,—,z, 0} while si^m € {+,—,z}. We intro¬ 
duce symbols “R” and “I” to distinguishes between real 
(0s-|-0|) and imaginary {iO^—iOD combinations of basis 
operators, respectively. We use also letters “E” and “O” 
to distinguishes between operators which respectively are 
odd and even under the spin-flip transformation. Since 
the Hamiltonian is invariant under spin-flip and time- 
reversal transformation, we separately study four orthog¬ 
onal sectors of operators denoted as RE (includes, e.g., 
the Hamiltonian), IE (includes, e.g., the energy current), 
10 (includes, e.g., the spin current), and RO. In the in- 
tegrable parent model, the local and quasilocal CQ exist 
in all four sectors provided A < Ii^ In order not to ex¬ 
clude any symmetry sector from our considerations, we 
take A = i. We restrict also the Hilbert space to the 
states with = 0. 

As follows from Eq, ®, breaking of integrability af¬ 
fects the correlation functions by either changing the sup¬ 
port of Qi (parameterized by A;) or by changing the pro¬ 
jected operators, Qf^ ■ In order to quantify the latter 
changes we have calculated the projections 


{Qfm 

wQnww 


( 12 ) 


where Qo,/ are CQ obtained for the integrable parent 
model. Figured] shows size-dependence of A; and Pi for 
various r. We observe that the perturbation strongly re¬ 
duces A/ (except for Ai in RE sector discussed below), 
whereas the projected operators do not change sig¬ 
nificantly. Therefore, we conclude that the main effect of 
perturbation consists in expanding the support of ACQ. 

From now on we focus on the support of ACQ as pa¬ 
rameterized by Aj [see Eq. ©]• The leading eigenvalues 
in the parent integrable model (lines with points in Fig. 
d]) are independent of the time-scale t indicating that 
the corresponding Qo,; are strictly conserved. It holds 
true both for local CQ with A/ = I (Figs. dJ^- anddb) as 
well as quasilocal CQ with 0 < A; < 1 (Fig. dt)- For 
the perturbed system, the only strictly conserved quan¬ 
tity is Qi in the RE sector, which actually represents 
the Hamiltonian. All other Qi are quasilocal for finite 
T. Their supports visibly depend on the time-scale even 
for quite large r, hence they are conserved only approxi¬ 
mately. We have verified these conclusions also for other 
eigenvalues, symmetry sectors, supports M = 3,4,5 and 
perturbations a/A = g, j. 

The most relevant and challenging problem is to estab¬ 
lish A; for large r and small a. It is also important that 
the finite size scaling L — >■ oo precedes the limit t —>■ ooi^ 
In Fig. [2] we plot A; linearly extrapolated to 1/L —>■ 0 
and normalized to results for integrable parent model. 
We clearly see that the dependence of extrapolated A/ on 
r and a is universal 


i?;(r,a) 


A/(L — >• oo, T, a) 

A;(L —>■ oo, r —>■ oo. Of = 0) 


Riira'^), 


(13) 



FIG. 1. Size-dependence of the leading eigenvalues A; [Eq. 
©1 and projections Pi [Eq. (I12ll 1 for various time-windows 
[0, r]. The symmetry sector (RE or 10) and other parameters 
are shown in each panel. Lines with points in panels a,b and 
c show results for the parent integrable model. Note different 
scale in d. 
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FIG. 2. Extrapolated and normalized eigenvalues A; [see Eq. 
1131 as a function of l/ro^. The vertical sizes of symbols 
show the differences between linear extrapolation from L = 
12,14,16,18 and L = 14,16,18. 


and can be well approximated by 

(t, a) ~ — arctan f—\— ). (14) 

TT / 

We have found the same type of behavior for other eigen¬ 
values (excluding the Hamiltonian), symmetry sectors 
and accessible supports M = 3,4,5 (not shown). 

Our main result for the Heisenberg model at infinite 
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temperature [Eqs. (IT3)) and (O] can be best explained 
in the formalism of the memory functions. We apply the 

S tion procedure according to Mori^ (see also Refs. 

, to analyze the relaxation function for CQ of the 
integrable parent model, Qq^i 


4>z(w) 

M{uj) 


(Qo,/l7r^lQo.0 = - 

L — uj w + M (u!) 

(P£Qo,il p^p_J P£Qo,i), 


(15) 


where LA = [H, A] and P is the projection onto the 
operator space orthogonal to Qo,i- 

When the formalism is applied to the perturbed inte¬ 
grable system (nni) with a <g; 1 , it follows directly from 
Eq. (fTSl) that 


LQo^i = a[H', Qo,i] oc a, (16) 


so that M{oj) = a^M.{w). It is plausible, but by no 
means obvious that the imaginary part of the memory 
function, is almost constant for small |w|. How¬ 

ever, if the latter is true then the dynamical relaxation 
reduces to 


4>z(w) ^ 


WQoA? 

UJ + ia'^M"{uj) 


IIQo.ziP 

UJ + ia^ji 


(17) 


We end up with a Lorentzian form which explains the 
specific dependence of eigenvalues A; on the time-scale r 
and perturbation a in Eq. (HI- Namely, using Eqs. 
and (El we find 

/ dw Im [$i(u;)] oc arctan ( —— ). (18) 

J-i 


The memory function analysis can be easily general¬ 
ized also to a quantity A which is not CQ but has sub¬ 
stantial overlap with the conserved quantity. 


U\Qo,)\ 

ll^ll IIQo./ 


(19) 


In such case the numerator in Eq. El) should be renor¬ 
malized becoming the Drude weight (dissipationless part) 
of the considered operator A 


_ ^ \\Qo,i\nA\Qo,i)\^ ^ \\Qo,\nA\Qo,0\^ 

^ UJia‘^(cu) a; + 

( 20 ) 

The latter equation is valid only for weak enough pertur¬ 
bation a I and in the low uj regime. 

For the numerical calculations of the memory function 
we employ the microcanonical Lanczos method,— well 
adapted for the studies of dynamics at /3 —>■ 0 where 
we can evaluate spin systems with up to L = 32 sites. 
The important parameter is the number of Lanczos steps 
Nl < 20000 which determines the u) resolution of the 
method 6uj ^ SE/N]^ where 6E is the energy span of the 
L-site spectrum, so that we reach 5uj ~ 10“^. 



FIG. 3. Panels a and b show normalized M" jo? for the en¬ 
ergy current jT, and M" for the spin current js, respectively. 
Points in panels c and d show the same as in Fig. 2 but for 
M = 3 (note that jE,ja £ Here, lines show left hand 

side of Eq. m for js (in c) and ja (in d). The latter quanti¬ 
ties have been normalized to the results for integrable parent 
model. 


We have carried out numerical calculations at /? —>■ 0 
for Qo,i in the IE sector being the energy current js in 
the unperturbed parent model and for the spin current 
js which has large projection on quasilocal Qo.i in the 
10 sector. From numerically obtained $"( 0 ;) we have 
extracted the relevant M{uj) via Eq. (ITSl) . Results pre¬ 
sented in Fig. confirm that M"{uj)/a^ for Je is in¬ 
deed very broad, featureless in a wide range uj € [ 0 ,a;o] 
uiQ > 1, and (almost) independent of a. On the other 
hand, $ 1 ( 0 ;) for js has a nonzero M"{uj) even for inte¬ 
grable a = 0 case, since js is not conserved. Still, in the 
regime w <C 1, M"[uj) as well <i);(w) follow the scaling as 
given by Eq. El- Note that for the Lorentzian assump¬ 
tion in Eq. El it is enough that it holds for w < 0 ^ 71 . 
Most importantly. Figs. [3]: and|3ji show convincing quan¬ 
titative agreement between the result obtained from our 
general approach and the formalism of the memory func¬ 
tions. Since the latter results have been obtained for 
much larger systems (but for two observables only) they 
can also serve a test of the finite-size scaling of A;. 


IV. DISCUSSION AND CONCLUSIONS 

We have proposed a general algorithm for a construc¬ 
tion and quantitative description of a full set of almost 
conserved quasi-local operators of weakly non-integrable 
lattice systems. The method has been implemented in 
the generic case of a XXZ model perturbed away from 
integrability by the 2 nd nearest-neighbor interaction of 
strength a. We have obtained a set of orthogonal ap- 
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proximately conserved operators, Qi. For an assumed 
time-scale r, these quantities completely determine the 
correlation functions of all local observables supported 
on several (M) lattice sites. Qi smoothly depend on r 
and a, and for the limiting case r —>■ oo, a —>■ 0 coin¬ 
cide with strictly conserved (local or quasilocal) quanti¬ 
ties of integrable parent model. We have shown that the 
perturbation influences the correlation functions mostly 
by expanding the supports of Qi. In our approach 
this effect is parameterized by eigenvalues A; decreas¬ 
ing from A; = 1 (for local operators) down to A; = 0 
(for generic nonlocal operators). We have found a scal¬ 
ing Ai/Af=° oc arctan [l/(ra^ 7 i)] universal for all Qi ex¬ 
cept the Hamiltonian. For a 7 ^ 0, the Hamiltonian re¬ 
mains the only strictly local or strictly conserved quantity 
(Ai = 1 independently of r) while all other Qi become 
quasilocal. We have found the same scaling also for a sys¬ 
tem perturbed by 4th nearest-neighbor interaction (not 
shown). This scaling seems to be typical, however, one 
cannot exclude that it breaks down for other specially 
tuned perturbations (see, e.g.,|43, and ls^ - ls^l . 

We have found a qualitative and quantitative agree¬ 
ment between our results and the memory functions ob¬ 
tained numerically for spin and energy currents for much 
larger systems. The latter analysis allowed us also to ex¬ 
plain the origin of the specific scaling of A/. Since this ex¬ 
planation is of perturbative character, we believe that the 
validity of the obtained scaling extends down to a —>■ 0 
well beyond the regime which can be inferred directly 
from bare numerical results. 

Within each symmetry sector we have found that the 
smaller A/ is (roughly understood as a more extended 
support of Qi) the larger is the scattering rate A 

relevant open question emerges: how many independent 
scattering rates are introduced by a single perturbation? 
Since all the scattering rates found in our studies are of 
the same order of magnitude, this problem may pose a 
challenge. 
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Appendix A: The choice of the scalar product 


As an alternative to the scalar product defined in Eq. 
© one may consider also other scalar products discussed. 


e.g. in Ref. [H 


{A\B), 


^(AtR + RAt) 

1 /t* D P{Bn)+PiEm) 

j -^mn-^rnn ^ ’ 

mn 


(Al) 


or 

(A|R )2 = ^ 

P-^ Jo 

= A*mnBmnP{En) >(^ 2 ) 

where A^n = (TO|A|n), H\n) = R„|n). Both these 
scalar products maintain the essential property, i.e., 
(A"|R")i = (A"|R)i and (A"|R ")2 = (A"|R) 2 . Cal¬ 
culating the scalar products of operators averaged over 
infinite time-window we find that the only contribution 
comes from states with equal energies (R„ = Rm), hence 

(A|R) = (A|R)i = (A|R)2 = ^ A*^^B^MEn). 

(A3) 

Consequently, the stiffness -^(AA) can be expressed 
in the same way by all considered scalar products 
as (A|A)(i 2 ). However, when discussing the memory 
function^ at finite temperature /3 < oo one should use 
the scalar product defined in Eq. (IA2I) . 


Appendix B: Support of the approximately 
conserved quantities 

The orthogonal set of ACQ consists of operators 


Qi=Y,UsiOl, (Bl) 

S 


where the unitary matrix U is defined in Eq. (jS]) and the 
norm of Qi can be found as 

WQiW" = J2(UsiOl\U,,idl) = J2uUoi\oi,)Usn = Xi. 

s,s' s,s' 

(B2) 

We split ACQ into two components Qi = + Q^, 

where the former operator is supported on M sites = 
Es'(Os'\Qi)Os' G Bff while the latter one {Q^\Os) = 0. 
Using Eq. (IBII) . Eq. ([5]) and the identity (A|R) = (A|R) 
we find 

= Y.iOs'\o:)UsiOs' = Y,{6:,\o:)UsiOs' 

s,s' s,s' 

= Y,XiUsnOs'. (B3) 
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The latter result together with the assumption concern¬ 
ing the orthonormal basis, (Os'l^s) = Sgs' yields 

IIQMP = (5] XiUs-iOs'\J2>^iUsiOs) 

s' s 

= XfY,ul,{Os-\Os)U,i = Xf (B4) 

s,s' 

Eqs. (jB2ll and (IB4p lead to Eq. 0 which relates eigen¬ 
value Xi with the support of Qi. 


Appendix C: Correlation functions and saturated 
Mazur bound 

We consider an operator supported on M sites, A G 
which can be expressed in terms of basis operators 


A = J2s Using Eq. (IBII) we find 

A^ = Y,ciso: = Y,^sUiQi = Y,yiQG (Cl) 


and 


{Qi\A) = {Qi\A^) = viXi. 


(C2) 


Repeating the same calculations for some other operator 
B G Bff we arrive at Eq. ([8|): 




= E 


l:Xi^0 


A/ 

{A\Qi){Qi\B) 

iQilQi) 


Xi 


(C3) 
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